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1. Introduction 

The Kadomtsev-Petviashvili equation in its version called KPI 

(lit - 6uu Xl + u XlXlX1 ) Xl = 3u X2X2 , (1) 

is a (2+l)-dimensional generalization of the celebrated Korteweg-de Vries (KdV) 
equation. As a consequence, the KPI equation admits solutions that behave at space 
infinity like the solutions of the KdV equation. For instance, if u\(t,Xi) obeys KdV, 
then u(t, x\, x%) = Ui(t, xi + ^ix% + 3fi 2 t) solves KPI for an arbitrary constant /z € R. 
Thus, it is natural to consider solutions of that are not decaying in all directions 
at space infinity but have 1-dimensional rays with behaviour of the type of u\. Even 
though KPI has been known to be integrable for about three decades 2 , :<| , its general 
theory is far from being complete. Indeed, the Cauchy problem for KPI with rapidly 
decaying initial data was resolved in by using the Inverse Scattering Transform 

(1ST) method on the base of the spectral analysis of the Nonstationary Schrodinger 
operator 

£(x,id x ) =id X2 +d 2 Xi -u(x), x=(x 1 ,x 2 ), (2) 

that gives the associated linear problem for the KPI equation. However, it is known 
that the standard approach to the spectral theory of the operator J2J), based on integral 
equations for the Jost solutions, fails for potentials with one-dimensional asymptotic 
behaviour. 

In ^OJ-JS] the method of the "extended resolvent" (or, for short, method of 
resolvent) was suggested as a way of pursuing a generalization of the 1ST that enables 
studying the spectral theory of operators with nontrivial asymptotic behaviour at 
space infinity. In f° r the Nonstationary Schrodinger and heat operators the 

case where there is only one direction of nondecaying behaviour was considered. The 
starting point in solving the problem was the embedding of the pure one-dimensional 
case in the two-dimensional spectral theory, building the two-dimensional extended 
resolvent for a potential u(x) = u±(xi). Then, a potential u(x) = Ui(xi) + v,2(x), 
where U2(x) is an arbitrary decaying smooth function of both spatial variables, was 
considered and the corresponding resolvent was constructed by dressing the above 
resolvent for u(x) — u%(xi). Finally, all mathematical entities generalizing the 
standard ones in 1ST, as Jost solutions and spectral data, were derived by a reduction 
procedure from this dressed resolvent. 

Here, we consider the case of a potential not decaying along multiple non parallel 
rays, which is substantially more complicated since we have not the one-dimensional 
sample as a guide to follow and we must construct directly the resolvent, without 
passing trough the embedding of one-dimensional entities in two dimensions. 

Therefore, we are obliged to consider directly true bidimensional potentials. In 
|2()| . by using recursively a binary Darboux transformation |21| . it was constructed 
explicitly not only a two dimensional potential u which describes N solitons [22] of 
the most general form "superimposed" to a generic background but also its Jost 
solutions. However, this recursive procedure seems not to be easily generalizable to 
the construction of the corresponding extended resolvent. On the other side, it is 
known (see [21]) that, in the framework of the extended resolvent approach, the whole 
hierarchy of time evolution equations related to £ in (J2J) , which can be considered as 
infinitesimal Darboux transformations, can be obtained by considering the similarity 
transformation L = QLC^ of the extended version L of £, where Q is a convenient 
unitary operator. Here, we show that, by using a twisted transformation L( — LC} 
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with C satisfying weaker conditions, one can bypass the recursive procedure and build 
directly the final potential u and Jost solutions. Then, we use this^operator that we 
call twisting operator, to build directly the extended resolvent of L as a bilinear form 
in terms of the Jost solutions. The main properties of the Jost and so called auxiliary 
Jost solutions and of the resolvent are studied. 

In a forthcoming paper, following the method developed in ^H] - CSj ; we generalize 
the result obtained in this paper. Precisely, we consider the potential obtained by 
adding to the potential describing N solitons an arbitrary bidimensional smooth 
perturbation, we construct the corresponding extended resolvent by dressing this one 
obtained in this paper and we derive the corresponding Inverse Scattering problem. 

2. Background theory 

2.1. Extended operators and resolvent 

In this section we briefly review the basic elements of the extended resolvent approach. 
For further details, we refer the interested readers to 

Let us consider the operators with kernel L(x,x') = C(x,id x )8(x — x') where 
jC(x,id x ) denotes a differential operator whose coefficients are smooth functions of x 
and let us introduce what we call the extension of these differential operators, i.e., 
to any differential operator C we associate the operator L(q) with kernel 

L(x, x'; q) = C(x, id x + q)S(x - x') = e 1 ^^ L{x, x'), (3) 

where x = (x\,X2), x' — {x'^x^) € M 2 , and q = (qi, q2) 6 C 2 and 

qa; = qi^i + q 2 £ 2 . 

The q variable will play in the following the role of a spectral parameter and we use 
a bold face character to emphasize that it is complex. By using the Fourier transform 
we can write 

L(x, x'\ q) = J dae- ta{x ~ x ">£{x,a + q), a=(a 1 ,a 2 ). 

Then, it is natural to introduce more general operators -A(q) with kernel 

A(x, x'; q) = J da e -*»C*-* V(x, a + q) (4) 

obtained by considering not just a polynomial C{x, q) in q but a tempered distribution 
P(x,q) of the six real variables x, qjf and qcj. Notice that 

A(x,x';c l ) = e i ^-^A(x,x';q), q = q a , (5) 

and that A(x,x';q) belong to the space S' of tempered distributions of the six real 
variables x, x' and q — (51,52)- Since definition up to the introduction of 
the spectral parameter q by shifting a, coincides with the definition of a pseudo- 
differential operator we shall call the operators belonging to the space S' extended 
pseudo-differential operators, or by short operators, and V(x, q) their symbol. 

In the following it results often useful to use instead of the symbol V(x, q) its 
Fourier transform with respect to x, i.e. 

^feq) = 7A2 / dxe ipx V(x,q), p= (pi,p 2 ). 
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From an d © it follows that A(j>; q) is related to A(x, x'; q) by 

A(p;q) = y dx y dxV (p+qsR)MC,3}a; 'A(x,x';q). (6) 

Then, we consider A(x ) x';q) and A(p;q) as the representation of the operator A(q), 
respectively, in the x-space and in the p-space. The inverse of (JSJ) is given by 

A(x, x'; q 9 ) = y dp y dq^ e -^+**)*+^'A(p; q). (7) 

On the space of this operators we define the hermitian conjugation as 

A\x,x';q)=A(x',x;-q), A\p; q) = A(-p;q + p), (8) 

in terms of kernels in x or p-spaces. For generic operators A(q) and B(q) with kernels 
A(x, x'; q) and B(x, x'; q) we introduce the standard composition law 

(AB)(x,x';q) = [ dx" A(x, x"; q) B(x", x'; q), (9) 



if the integral exists in terms of distributions. In terms of kernels A(p; q) and B(p; q) 
this composition takes the form of a shifted convolution 

(AB)(p;q) = jdp'Aip-p'^ + p^Blp'-^). (10) 

An operator A can have an inverse A -1 in the sense of this composition, i.e., such 
that AA~ X = I or A^ 1 A = I, where I is the unity operator in S 1 , I(x, x'; q) = S(x — x'), 
being S(x) = 5(x\)S(x2) the two-dimensional 5-function (or I(p; q) = S(p) in p-space). 

Of course the two representations in the x and in the p-space are equivalent and, 
in principle, one could work always in one of them. However, it results often convenient 
to pass from one representation to the other. Thus, the p-space is more suitable to 
study analyticity properties, while boundedness is more easily studied in the x-space. 

The extension of the Nonstationary Schrodinger operator J2J is given by 

L = L Q -U (11) 
where in the x-space 

L Q (x, x'\ q) = [i(d X2 + q 2 ) + (d Xl + qi ) 2 ] 6(x - x'), U(x, x'; q) = u(x)5(x - x'). (12) 

The main object of our approach is the extended resolvent (or resolvent for short) 
M(q) of the operator L(q), which is defined as the inverse of the operator L, i.e., 

LM = ML = I. (13) 

Here we omit to specify the additional conditions that guarantee uniqueness of the 
extended resolvent as solution of H13[) . referring, say, to PIUS] . 

For a real potential u(x), as we always consider in the following, we have 

= L, Af f = M. (14) 

Now, let us consider in this section the case of a rapidly decaying potential u(x) 
in (P)l. 

One of the main advantage of the resolvent approach is that the dressing operators 
can be obtained directly from the resolvent by means of a truncation and reduction 
procedure. Thus, here the dressing operators v and lu are defined by 

i/(p;qi) = (MLo)(p;q) , , w(p;qi) = (L M){p;q) , (15) 



e(qi) 



(qi+pi)-p 
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where we introduced the special two-component vector 

£(k) = (k,k 2 ). (16) 

In fact these operators dress, in the sense proposed and developed by Zakharov- 
Shabat 0, the operator L(q) and its resolvent M(q). More precisely, they admit the 
following bilinear representation in terms of v and to 

L = vL uj, M = vM w, (17) 

where L and M are the bare operators 

L (p;q) = 6(p)(q 2 -q 2 1 ), M (p;q) = -^-j. (18) 

q2 - qi 

Dressing operators are mutually adjoint 

v i = u, (19) 

mutually inverse 

wv = I, (20) 
vu = I, (21) 
and obey the equations 

Lv = vLq, ujL = Lquj. (22) 
The kernels of these operators in p-space obey asymptotic 

lim v(p;q) = S(p), lim ui(p;q) = S(p), (23) 

qi — too qi — ^oo 

are independent of q 2 and analytic functions of the variable qi in the upper and lower 
half planes. 

Spectral data and the inverse problem can be formulated in this operatorial 
approach. Here, we give not details but only the main formulae, that will be useful 
in the following. Let v^- and uj^ denote the operators with kernels being the limiting 
values of v(p\ q) and u>(p; q) at the real qi-axis from above and below 

v ± {p;c\x)=v{jp;q m ±iQ), w ± (p; q : ) = w(p; qi» ± iO). (24) 
Then we get the relations 

v ± = v T F T 1 lu ± =F ± uj t , (25) 
where we introduced the spectral data 

F ± =lo ± V t. (26) 

By construction the kernels F^ (p; q) depend on three real variables p\ , pi , and qigj 
and thanks to (|19() - (|22|l these operators obey 

L (q)F ± = F ± L Q {q), when Ql = 0, (27) 
(^ ± )t=F ± , (28) 
F+F- = I. (29) 

By means of the last equality and (|23|l it can be shown that the kernels F^ 1 (p; q) have 
the representation 

F ± {p;<l) = ${p) +${P2 ~Pi(pi + 2qisff))/ ± (pi +qiSR,qisre), (30) 
where the functions / ± (pi,qisR) depend only on two real variables. 
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2.2. Hat-kernels, Green's functions and Jost solutions 

Let us associate to any operator A(q) with kernel A(x,x';q) its "hat-kernel" 

A(x, x';q) = e^'^ A(x, x';q). (31) 

If A{q) is an extended differential operator L(q) then this procedure is the inverse 
of i.e., L(x,x';q) — C(x,x'), while for a generic operator the hat-kernel can 
continue to depend on q. For any extended differential operator L(q) and for any (not 
necessary differential) operator B(q) the following relations hold 

(LB)(x, x';q) = C(x, d x )B(x, x'; q), 

(BL)(x,x';q) = C d {x\ d x ,)B(x, x'; q), 

where C A is the operator dual to C In particular, by 1131) we have that 

C(x,d x )M(x,x';q) = C d (x', d x ,)M(x, x'; q) = S(x - x'), 

so that the hat-kernel of the resolvent defines a family of Green's function depending 
on the bidimensional parameter q. In particular, the Jost solutions are defined by 
means of the following Green's function depending on a complex spectral parameter 
k 

g(x : x',-k) = M(x,x';lcs(k)), (32) 

with the vector £(k) defined as in (|16l) . while the advanced/retarded solutions are 
defined by the Green's functions Q±{x, x') obtained by the following limiting procedure 

G±(x,x')— lim lim M(x,x';q). (33) 
The Jost solution $(x,k) and its dual ^(a;,k) are defined by 

k) = J dx'e- u{u)x 'c d {x\ d x ,)G(x, x\ k) (34) 

<&(x',k) = J dxe u ^ x C {x,d x )g(x,x',k), (35) 

or by using definitions l|15[) in terms of the dressing operators by 

$(x, k) = e- u ^ )x X (x, k), V(x, k) = e u W x £(x, k), (36) 

where 

X (a:,qi) - J dpe^u^qx), £(z, qi ) = J dpe-^wfam - Pl ). (37) 

The q2 independence of the kernels v{p\ q) and u>(p; q) implies, thanks to Q and l|31|) . 
that the corresponding hat-kernels v{x, x'; q) and u>(x, x'; q) are independent of q2 and 
are proportional to S(x2 — x' 2 ). Thus, instead of 1)36(1 and (|37|l we get relations 

$(z, qi ) = J dyV(x,y;q^)e' a ^y, tffoqi) = J dye^^y^;^), (38) 

Writing (T3J in the form M = M + M UM one can derive from l|34|) and (|35|l 

the standard integral equation for the Jost solution in the case of a smooth potential 
rapidly decaying at space infinity 



(39) 
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where 

G (x,x',k) = ^-J da9(a^ X2 )e^^\ 

is the Green's function, introduced in [5]. In our approach this Green's function follows 
directly from the second equality in Ijl8|) . transformation J2J, and definition l|32(l . 
Solvability of this integral equation under some small norm assumptions was proved 
in jB] and thanks to 139J1 it is easy to show that x( x , k) has the asymptotic behaviour 

lim xO,k) = 1. (40) 

x — >oo 

Properties of the dressing operators lead to corresponding properties of the Jost 
solutions. Thus flUI gives 



$(z,k) = tf(z,k), X (z,k)=£(z,k), keC. (41) 
Property is equivalent to 

lim x(;r,k) = 1, lim £{x,k) = 1, (42) 

k — >oc k — >oo 

relations (|2*U|) and (|2*l"1) can be considered as scalar product and completeness of the 
Jost solutions: 



n — / dxi$(x,k)*(x,k + a) = 5(a), keC, a € R, (43) 

-!- / dk K $(z, k)tf (j/, k) = S( Xl - yi ); (44) 

and equalities l|22|l are equivalent to the Nonstationary Schrodinger equation and its 
dual: 

C(x,d x )$(x,k) = 0, C d (x,d x )V(x,k)=0. (45) 

The Jost solutions are analytic in the complex plane of the spectral parameter k, 
ka 0. Using for the boundary values at the real k-axis notations of the type (|24|l 
we get from (I25|l and i|30|) the standard 5 nonlocal Riemann-Hilbert problem 

$ ± (a;,fc) = $ T (x,k) +J da<f> T (x,a)f T (a,k), keR. (46) 



3. Twisting transformation 



Here, we construct a two dimensional potential together with its Jost solutions which 
describes N solitons |22l 12^] "superimposed" to a generic background by using a 
twisted transformation L( = QL from the extended differential operator L in 
to a new operator L of the same form obtained by means of an operator Q which 
is isometric but not necessarily unitary. The use of the twisting operator Q allows 
us, bypassing the usual procedure consisting in applying recursively binary Darboux 
transformations, to get directly not only the Jost and auxiliary Jost solutions but also 
the extended resolvent, which results to be a bilinear expression in terms of these Jost 
and auxiliary Jost solutions. 



On the extended resolvent of the Nonstationary Schrddinger operator 8 
3.1. Properties of twisting operator C, 

Let us consider twisting the operator L in to a new operator L of the same kind 
L = L -U, U{x,x':q) = u{x)5{x - x'), (47) 

by means of an operator £ according to the formula 

L( = (L. (48) 

We consider a potential u(x) in L which is real, smooth and rapidly decaying at space 
infinity and we search for a £ such that the new potential u(x) is also real and smooth, 
while condition of rapid decaying is not imposed. Notice that since L and L are both 
selfadjoint from l|48l) it follows that 

<?L = L<?. (49) 
In addition we require that £ obeys the conditions 

(I) the operator £ is isometric 

C f C = I, (50) 
but not necessarily unitary; 
(II) the kernel £(p; q) is independent of q2, 

C(p;q) = C(?;qx); (si) 

(III) the kernel ((p; q) obeys the asymptotic condition 

lim C(p;q)=«(p)- (52) 

qi— >oo 

Then, a specific transformation £ is chosen by fixing its analyticity properties in 

qi- 

Notice that because C is not unitary the operator 

P = /-CC f , (53) 
is not zero and since it satisfies 

=P (54) 

and, thanks to (|5Up . 

P 2 = P, (55) 

it is an orthogonal projector. 

We fix the analyticity properties of £ by requiring that it generates not only the 
new potential u via (|48|l but also the new dressing operators v and w. In fact, taking 
into account we get by igSJ and i|4"9")) 

LC,v = C,vL , u(?L = L Q utf, 
and, therefore, the operators z/ and w defined by 

i7r = Ci/, (56) 

r^^wC 1 , (57) 
for any operator r commuting with Lo satisfy the equations 

Lv = vLq, ojL — Lquj. (58) 
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analogous to the equations lll'l'l) satisfied by dressing operators v and ui of L. Since v 
and u> are are mutually adjoint 

i? + =2, (59) 

we can limit ourself to consider v and say that v can be considered the new dressing 
operator of L if the operator r and consequently Q is chosen in such a way that the 
kernel of the v in p-space 

(IV) independent of q 2 

£(p;q) = v(p;qi), (60) 

(V) is an analytic function of the variable qi in the upper and lower half planes, 
continuous on the two sides of the real axis and obeys the asymptotic behaviour 

lim v{p] q) = 5(p). (61) 

Since v, u>, and £ by ijTIll . are independent of q 2 we deduce from (IIVII that also r 
must be independent of q 2 and, then, since r commutes with L , that its kernel has 
the form 

r(p;q) = *(p)T(qi). (62) 
Moreover, from l|V|l we get that r(qi) has asymptotic behaviour 

lim r(qi) = 1. (63) 

Taking into account l|2l)|l and (|50|l we get from l|56|l and (|57H the scalar product of the 
new dressing operators 

wv = T'\ (64) 

where we introduced the selfadjoint operator 

T = tt^ (65) 

with kernel 

T(p;q)=*(p)*(qi) (66) 
where thanks to the composition law (|10fl and definition (|SJ 



*(qi)=T(q 1 )r(q 1 ). (67) 

From l|64|) . thanks to JVJ, we deduce that l/t(qi) and therefore l/r(qi) must be 
chosen to be a function analytic in the upper and lower half qi-planes and continuous 
on the two sides of the real axis. Therefore, we assume that 

(VI) r(qi) is meromorphic in the upper and lower half planes without zeros and with 
a finite number of poles and satisfies asymptotic (|63|) . 

Consequently, t(qi), which plays the role of transmission coefficient, is analytic 
in the upper and lower half qi-planes, continuous on the two sides of the real axis 
with no zeros and with poles at the poles of r(qi) and r^). 

The completeness relation of the new dressing operators, which will play a crucial 
role in the following, can be written in terms of T and the projection operator P. In 
fact from JSSJl, (JSJ||, and JHSJl we get 

VTZj + P = I. (68) 
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It is also worth to mention that the operator P annihilates the new dressing operators. 
Indeed, from Q64|l and l|68|l we have 

Pv = ZjP = 0. (69) 

Finally, by using the completeness relation l|21(l from Ij56|l we have 

C = vtu (70) 

and, once given an operator r satisfying the above described conditions, the analyticity 
properties of C, are fixed and consequently £ itself. However, in determining the 
analyticity properties of £ we have not yet required that also the new dressing operators 
v and ui satisfy non local Riemann-Hilbert problems analogous to them satisfied by v 
and to in (|25|l . This will be done in the next section. 



3.2. Transformed continuous spectrum 

In order to describe the discontinuity at the real axis in the complex qi plane of 
the twisting operator £ and of the new dressing operators and Jost solutions we use 
notations (|24[1 and mention that for any operator A by definition © we have 

(A^ = (A^. (71) 
From i|56|) m the limits qi — » ±0 we get 

u ± r ± = ( ± v ± . (72) 
Multiplying from the left by (C ) , thanks to Ij50|l . we get 

u ± = (C t ) ± ? ± T ± (73) 

Inserting H25fl into the r.h.s. of l|72|l and then inserting v T from (|73|l in the obtained 
equation we get 

IrM = C ± (C t ) T ^ T ^ T ^ T - 

If we impose that v satisfies a Riemann-Hilbert problem analogous to that one in (|25|l 
satisfied by v we deduce, taking into account that according to I|tj9|) v T has a left 
annihilator P^ , that 

C ± (C t ) T = Cl7 + C 2 (/ + P=F) = ClI + C2C T (C t)T 

with ci and c 2 constant. Then, since C satisfies the asymptotic property (|52|) thanks 
to l(5rJ|) we derive 

C + =C, (74) 

that is 

(VII) the twisting operator £ is continuous on the real axis of the qi plane. 

Let us, now, introduce the operator 

p± =T ± J? ±( r ±)t. (75) 

Its kernel equals (see (|20J>) 

F^faq) = «5(p)|T ± (qi 8f i)| 2 + 6(p 2 -Pl{pi + 2q 1 s R ))/ ± (p 1 +qi5R,q m ), (76) 
where 

7 ± (pi,qi») =T ± (p 1 )/ ± (p 1 ,q lffi )7±(q 15R ). (77) 
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Thanks to properties of r , the operators have properties l|27|) and (|28[1 . Taking 
into account (|71|l and l|(jtj|) we get finally 

v^T* = v T F T , T ± uj ± = F^uF, (78) 
where the second equation is obtained by conjugation. Comparison of l|78|l with (|25fl 
motivates the definition of F^ as the continuous spectrum of the new potential u. 
Again, thanks to (|71l) and l|6(jf> . we see that relations (|29[) and <|2(i|) are modified, 
respectively, as 

(T ± )- 1 J F ± (T =F )- 1 F =F = /, (79) 

and as 

= (T ± )- 1 J F ± (T^)- 1 . (80) 

In the case where the background potential u(x) is identically equal to zero the 
second term in (|3(Jll is absent and F^ = I. It is natural to preserve this property for 
the new spectral data, so we impose 

(VIII) unitarity: 

Ir^qm)! = 1, (81) 



that in terms of operators means that 
(r ±)t = ( T ±)-\ 

In analogy with 136|) and (|37|l we define 
X(x, qi ) = J dpe-^ufaqx), f(x, qi ) = J dpe^u^cn - Pl ), (82) 

k) = e- a ^ x x(x, k), *(as, k) = e a <»*Z(x, k), (83) 

where notation (|16|l is used. These functions can be defined also by relations 
generalizing (|38|l . Indeed, taking into account that the function r(qi) has no 
zeroes we have by ||5d]) and (|57|) that z/(p;qi) = (C^)(p; q)/ T (qi) and u;(p;qi) = 
(w(^)(p; q)/ T (pi + q~i)- Inserting this equalities in we derive by ©, lEU, an d (|3"S| 
that 

$(a;,qi) = — ^— / C(x, y; qio)$(zj, qi), (84) 

T (qi) J 

*(x,qi) = == / dyy(y,qi)(i(y,x;ciiQ). (85) 
T(qi) 

In what follows we prove that $ is the Jost solution of operator C and ^ is the 
Jost solution of the dual operator. Thanks to Ij59(l we have like in (|41|) 



$(a:,k) = *(x,k), x(x,k) = f(a:,k), k e C. (86) 

Under condition (Vj functions $(x,k), ^(x,k), ^(a;, k), and £(x, k) are analytic with 
respect to the variable k G C with only discontinuity at the real axis and functions 
x(x, k), £(x,k) obey asymptotic condition (021 • Eq. i(78|l enables us to write down 
relations between boundary values ^(x,k) and 4 ,± (x,fc) of these functions on the 
real axis. Indeed, taking into account definition ifTUI) of composition we get by ||BSJ|, 
(J7HJ), and JHU 

^(x,^^) = $T(x,k)+ J da$T{x,a)f T {a,k), k £ R, (87) 
that modifies relation 146|l . valid in the case of rapidly decaying potential u{x). 
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3. 3. Twisting transformation of the resolvent 

The operator C, once obtained the transformed operator L, can also be used to get 
the corresponding resolvent M = L^ 1 . From 148|) and the definition of the resolvent 
in (|13|l we have that M and M are related by the same relation as L and L, i.e., 



M( = CM. 

Multiplying l|48|) and 188|) from the right by ^ and recalling the definition (|53[) of P 
we get 

L = (Ltf+L A , M S MC • ,\/ A (89) 

with 

L A = LP = PL, M A = MP = PM, (90) 

where the second terms in both equations follow by hermitian conjugation. Since P 
is a projection operator we get directly 

M A P = PM A = M A , (91) 
L A M A = M A L A = P. (92) 

Viceversa, for Ma satisfying (|91|l and (|92|l the operator M in (|89|l is the resolvent of 
Z. In fact we have by (JHOJ), lO, and 

Zm = (CLC f + i A )(CMC t + M A ) =1 + CLC}M A + L A (M(K (93) 

To prove that the last two terms are zero we, first, note that thanks to lf5U|l and (1531 
we have that 

PC = ?p = o, 

and, then, that by ^ ( f M A = ( f PM = and as well L A C = LP( = 0. So 

ZM = ML = I, (94) 

where the second equation is obtained by conjugation. 

In conclusion, in order to obtain Ma, and consequently M via l|89|) . one can, first, 
find the general solution X (belonging to our space of operators) of the system 

PX = X, XP = X, (95) 

and, then, Ma will be the special X satisfying 

L A X = XL A = P. (96) 

3.4- Construction of operator C 

In correspondence with condition 1VI|) let r(qi) have poles at qi = Ai, . . . , Ajv and 
let us suppose, for simplicity, that they are simple and that 

A m 7^ A rl , A m ^ An, Vm/n 

A„3^0, Vn. (97) 
From properties (|VI|I and (jVIIIII of r by dispersion relation we get that 

r(qi ) = TTf qi ~-^ ^TT qi " — + z|A _" a|sgnqia , (98) 

„_! \qi - A„ / ^J; qi - A„ 
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For the residua of this function we have 



N /T . x e(A m3 A„ Q ) 

Am — A r , 



TV, 



res r( qi ) = -2iA m9 JJ m " . (99) 



1 7i=l, n7=m x 



From we get for the transmission coefficient 

^ \ qi - A„ / V qi - A n3 } - * I An»| / 

The residua at the poles at qi = A m and qi = A m are given by 

t m = res t(q 1 ) = 2zA m9 IT ( 7" ~ 7 ) , (101) 

qi=A m \A m ~ A n / 

71=1, n^m 

res i(qi) = t m . (102) 
qi=A m 

An alternative expression is given by 

t m = r(X m )T m . (103) 

Now, we have all we need for building the twisting operator £. First, we use (|2*TJl 
to rewrite (|5f>|t as 

C = vtuj. (104) 

Now, from the analyticity properties in qi of v stated in pVI) . ifVjl . of r as given 
in (|98fl and of w, from the continuity of £ on the real qi-axis stated in IjVIIf) and 

the asymptotic property of ((p; qi) we have that ((p; q) satisfies the following integral 
representation 

> / \ n , f , ,v(p-p';\ n )uj{p';X n -p[) 

C(p; q) = <5(p) + > r„ / dp = . (105) 

The kernel £(p; q) results to be piecewise analytic with discontinuities along the lines 
qiQ = — A n Q and it results do depend only on the values of ?(p;qi) at qi = A„. 
Therefore, in order to get £ it is sufficient to construct the new dressing operator v at 
the special values of poles of r. 

The kernel of C in the x-space is given by (0- Then, for its hat-kernel (see (|31Jl ) 
we get 

£(x, x'; q) = S(x — x') — i sgn(xi — x' 1 )5{x2 — x' 2 ) 

N 

xY / T nO{(qi+X n Q){x 1 -x' 1 ))^{x,Xn) 1 f{x',Xn), (106) 
n=l 

where we used notations (J2IIJ), (|3T|l . (fH^f) . and (|S3)l . By conjugation (see ©, lETTl) . 
and 186I) - ) we derive 

C^(x, a;'; g) = S(x — x') — i sgn(xi — x' 1 )5(x2 — x' 2 ) 

N 

x^r„0((qi - A^)^ -x' 1 ))$(x,A„)$(x',A„). (107) 

n=l 
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In order to complete the construction of £ we have to impose the isometry 
condition (|5(ip. Using i|10ti[l and I|1U7|) we get that this condition is equivalent to 



N 



S(x 2 - x' 2 )^f n e((qi - A n3 )(xi -xi))$(x,A n )*(x',A n ) 

71=1 

AT 

+ <5(x 2 - a;' 2 ) ^ T n 6((qx + \ n %)(xi - x[))$(x, X n )^(x', A„) 

n=l 

iV 

= i sgn(xi - a;i)5(x 2 - x 2 ) ^ T m T n ®(x, A m )*(x', A„) sgn(qi - A mS ) sgn(gi + A„ a ) 

m,n— 1 

x / dyi6»((gi - A m9 )(xi - yi))0((gi + A nQ )(yi - x' 1 ))$(y, A m )$(y, A n ) 



Let us, now, introduce the matrix 

(A m +A„) a oo 

0(i) 



dyi*(2/,A m )$(y,A n ) 



V2=X2 



N 



, (108) 



where the factor (A m + A„)q in the limit of integration defines the sign of the infinity. 
Thanks to (|83|) it is easy to check that these integrals are well defined if £(x, A m ) and 
x{x, A n ) are bounded at space infinity, and that by Ij86|l this matrix is hermitian: 

e m „(a:) = e mn (x). (109) 

Thus 

T m T n ^>(x, A m )$(x, A n ) = -9 a;i e mn (a;), (110) 
and the above condition is simplified to 

N 

^2 e (( a l - Am3)(.Tl - Am) 



/V 



T m *(a;', A m ) - i ^ 8 m „(a;')*(a;', A„) 



JV 



^ #((<7i + A ma )(a; 1 - a^)) 



JV 



T m $(x, A m ) + i ^ ®( X > ^n)Onm(x) 



V(x',\ m ) = 0, 



where a; 2 = x' 2 . The function in the l.h.s. is piecewise constant with respect to qi, so 
this condition is equivalent to the equalities 

~ 1 N 

<f>(x,\ m ) = - — V \„)Q nm (x), (111) 



n=l 
JV 



*(x, A TO ) = — ^ mn (x)*(x, A„). 



(112) 



By <|106[) in order to construct the hat kernel of the operator £ we have to 
determine the functions $(A m ) and ^(A m ). Thanks to (|lllf> and 11121) this means 
that we have to determine the matrix mn - For this sake we insert these equations 
in (TIOll , that gives 



^mn (^) ^ ^ ®mm' (x)^& (x, A m ' ) < l > (x, A n ' )0n'n (x) , 



On the extended resolvent of the Nonstationary Schrddinger operator 15 

that under assumption of invertibility of the matrix <d(x) can be rewritten as 

d Xl Q^ n {x) = *{x,X m )$(x, A n ). (113) 
Let us introduce in analogy with (|108[l the matrix 



B mn {x)= / dyi^(y,X m )^>(y,X n ) , (114) 

J V2=X-2 

-(A m +A„) a oo 

where again the limits of integration are uniquely determined by the asymptotic 
behaviour of A TO ) and *ff(x, X n ) given by ffity and ijlHJl. Moreover, 

-i 



lim B mn (x)e 



-i(e(x m )-e(x n ))x _ 



Xrn Xn 

and by (|41|l this matrix is hermitian, 



^mn(^) — -^nm(«E)- 

Now l|113fl gives 

Q mn (x) = (B(x)+C)^ n , (115) 

where we introduced the matrix C — ||c mn ||^ n=1 with matrix elements independent 
of X\. 

Let C± denote the matrices constructed from C as 

C± = Hew; m,n = l,...,N; ±A mS > 0, ±A„ Q > 0]|, (116) 

and C+ = / (C- = I) if all A„s are negative (positive). In |20I it was shown that if 
the matrix C is hermitian and obeys the positiveness condition 

±C±>0 (117) 

the determinant det(-B(a;) + C) has no zeros on the x-plane including infinity (for 
the case of zero background potential this result was obtained in |23|). Thus, the 
matrix Q(x) given by 1115fl exists for any x and the r.h.s.'s of 1111|) and 11121) are 
given explicitly in terms of the Jost solutions of the background potential u(x). Such 
matrix <d(x) is hermitian indeed, so that by (|41|l equalities l|lll|) and (|112l) give a 
special case of (I5BT) for k = X m . 

Finally, inserting 11111) in H06|) and l|112l) in (|107|l we get the hat-kernels of £, 
£t. Multiplying them by e~ q ( x ~ x ) (see (j2U) we get kernels 

((x, x'; q) = 5(x - x') - sgn^! - x\)S(x 2 - x' 2 )e- q ^ Xl ~ x '^ 

N 

x 6l (('?i + -aj' 1 ))*(a;,A m )e mn (a;)*(a; / ,A n ), (118) 



(1(x,x';q) = S{x - x') + sgn( Xl - x' 1 )S{x 2 - x' 2 )e- qi( - Xl - x '^ 

N 

x d (( qi ~ X m&)(xi -x[))<S>{x,X m )<d mn (x / )t>(x',X n ) 7 (119) 

m,n— 1 

proving by l|41(l and j 109(1 that they are mutually conjugate in the sense of (JHl- 
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3.5. Transformed dressing operators and Jost solutions 



By using the constructed twisting operator we can obtain the dressing operators 
by H5(j|) and l|57l) . while expressions for functions introduced in 1(83(1 follow from (|84|l 
and 1(85(1 . Indeed, inserting l|118|) and 1119(1 in these relations, we get 



$(a:,k) 

$(x,k) 
where we denoted 

so that by (j!14jl 



1 



r(k) 
1 



<J>(a;,k) 



-00 L 



*(x,k) - 



JV 

E 

m,n— 1 
N 

E 

m,n— 1 



A m )e m „(cc)/3 n (a;, k) 



/3 m (x, k)6 mn (a;)*(a;, A„) 



dj/i*(j/,A n )$(y,k) 



(120) 



(121) 



(122) 



V2=X2 



-(k+A„) a oo 



-^mn(x) — /^m(x, A n ) — /3 n (x, A m ). 

Thanks to the properties of the Jost solutions $ and \&, the functions /3 n (k) are 
analytic in the complex domain of k with exception of the real axis and of a pole at 
k = A„ , where by 

res /3„(x, k) = i. 

k=A„ 

By this proves that <I>(x, A m ) ($(x, A m )) given in 1111(1 (correspondingly, ((112(1 ) 
are values at k = A m (at k = A m ) of ((120(1 (correspondingly, 1(121(1 ). 

Taking into account the well known property of the determinants of bordered 
matrices 



1 



det r. 



r, 



7n+l,n+l 



r r _1 r 

ro+l,* n *,n+l> 



,n+l 

n+1,* Tn+l,n+l 

it is easy to see that ((120(1 is exactly the Jost solution constructed in [201 as result of 
N successive "binary" Backlund transformations. It is annulated by the operator 



C(x,id x )=id X2 +di 1 -u(x), (123) 

iff the matrix C in 1(115(1 is independent also on x-i. Correspondingly, ^(x,k) is the 
Jost solution of the dual operator. In 1(123(1 the new potential (cf. (J2J) is given by 
means of the well known [2] relation 

u(x) =u(a;) + 2<92 i logdete(izi). (124) 

We proved in [23j that this potential is smooth, real and finite for all x under 
condition 1(117(1 and that it decays in all directions on the x-plane with exception of 
a finite number of directions x\ — 2Ajsrie2 =const, where it tends to a unidimensional 
soliton. It was also proved there that x(x,k) (see (183(1 ) is a bounded function of its 
variables, analytic in k with a discontinuity at the real axis, and obeying asymptotic 
condition ((42(1 . while instead of ((40(1 we have 



lim 



xO,k) = 1, 



lim 



X0,k) = — . 

H k ) 



(125) 
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The Jost solutions given by (|120fl and i|121|) obviously obey l|8t)[) . Then properties of 
£(x, k) are the same up to the asymptotic 

lim f(x,k) = -l-, lim £(x,k) = l, (126) 

where we used that by (|u7|) t(k) — t(k). 

Summarizing, we have that kernels £(x,x';g), ^ (x,x';q) and v{p;c\), 5(p;q) (as 
given by (JH2Jl) belong to the space 5'(R 6 ), so they define operators in the sense of 
definitions given in Sec. [21 Moreover, these operators obey all conditions (|H)l ifVl 
that were imposed. Asymptotic behaviour (|125(l and (|126|l shows that the values of 
the Jost solutions at the poles of i(k), like in the one-dimensional case, have special 
relevance. Thus in addition to the values given in 1|) we have to consider <£>(x, Aj). 
Thanks to l(T2^|) and !jTT5|) we get from (fT^Ufl 

1 N 

$(x,Aj)=— — r ^ ®( X ^m)Q m n(x)c mj . 
3 m,n=l 

then by ifM^) . (fTH^I . and fTTTf> . ITT^I 

. N 

$(x, A m ) = — ^ ®( x > y^Cnm^, (127) 
tm n=l 

~ 1 N 

*(x, A m ) = r m c m „r n *(x, A,,). (128) 

n=l 

Using these equalities we get 

N 

^{t n $(x,A„)*(x',A„) + I n $(x, A„)*(x',A„)} = 0. (129) 

n=l 

Relation (|127f) together with (|87|) and the first equality in (|42|l close the formulation 
of the inverse problem for the Jost solution <fr(x, k). Analogously, l|128|) . l|57|) and the 
second equality in (|42|l give the inverse problem for the dual Jost solution ^(x, k). 

Boundedness o/<I>(x, A„), $(x, A„) and £/ieir dwa/ 

Here we prove that the Jost solutions <J>(x, A„), $(x, A„) and \t(x, A„), ^(x, A„) are 
bounded when x tends to infinity, while specific asymptotic behaviour essentially 
depends on the direction of the limiting procedure. Let us notice that ^(x, A„) and 
$(x, A„) can be written as 

$(x, A„) = e M^)*-;MA„)*~ (a;) Xn ), $(x, A„) = e- e *^+ u ^ x x(x,K), 

where x(x, A„) and x(x, A„) are known to be bounded. In the limit x — > oo along 
a direction of the x-plane some of the exponents e ^^ n ' x are increasing, some are 
decreasing or bounded. Taking into account that new potential, Jost solutions, and 
spectral data by construction are symmetric functions of Ai, . . . , Xn, we renumber 
them in a way that, say e iss ^ n ' x for n — 1, . . . , s — 1 are decreasing or constant and 
e ^(\n)x f or n, = S) . . . ( _/V are increasing, where s = 1, . . . , iV+1 is a number depending 
on the direction on the x-plane. Then $(x, A„) are decreasing for n = s, . . . , N and 
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A„) are bounded when n = 1, . . . , s — 1. In order to consider the complementary 
intervals we write H127J) for m > s in the form 

N s-1 

^ $(x, A„)r„c„ m = A m ) - ^ ®( x i ^n) r «c nm - 

n—s n— 1 

Thanks to conditions (|117fl matrix C cannot have zero eigenvalue. Then the same is 
valid for the matrix ||c mn ||^ n=s in the l.h.s. of this equality. Taking into account that 

all terms in the r.h.s. are bounded, we conclude that <fr(x, X n ) are bounded also in the 
interval n = s, . . . , N and then for all n. Now by (|127(l the same is valid for A n ). 
Boundedness of ^(x, A„), ^(x, A„) follows by conjugation. 



4. Transformed resolvent 



4.1. Operator P and completeness relation 



We already mentioned in Sec. VS. li the essential role played by the operator P. Here we 
derive an explicit expression for its kernel and present its properties. Inserting UlUtjfl 
and I1U7|) in (JHSJ and using (|114[) we get 



N 



P{x,x'\q) =isgn(x 1 - x' 1 )5(x 2 - x' 2 ) ^ T m#((<7i +A m s)(x 1 - x[)) 



x$(x, A m ) 



JV 



f(i',A m ) + i^T„$(i',A n )B mn (a;') 



n=l 
N 



+ isgn(x 1 - x' l )S(x 2 - x' 2 ) ^ T nO((qi - A„s)(a; 1 - x[)) 



71=1 



N 



$(x,\ n ) 



*(x',\ n ) 



By JTTH) A„) = i YZ=i T m®{x, A m )( 6(x))-^, so thanks to JEEJ and (JMl 



A„) - i r m $(x, \ m )B mn (x) = rr-<J>(a;, A„ 

m=l 

Thanks to this equality and its complex conjugate we get 
P(x, x'\ q) — iS(x 2 — x'2) sgn(x 1 — x[) 

N 

x [tnOdQi - A„s$)(a; 1 - x' x ))^{x, A„)*(x', A„) 

n=l 

+ ?n6 , ((9i + Kq){x 1 - xi))$(x, A„)*(x', A„) 

The form of this expression shows a discontinuity of the r.h.s. at x x — x[, while 
thanks to (|129fl the actual discontinuity is absent. In order to exploit this fact directly 
we write 

P(x, x'; q) = iS(x 2 - x' 2 )6(qx) 

N 

x [6{\n%)0{\(li\ - lA^D^^i - x[) + 9(-X n ^)9(x 1 - x[) 

n=l 
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- 0(A nO )0(|A n c*| - \qi\)e(x[ - Xl j\ $(x, X n )V(x', X n ) 
+ t n 6(\ n cs)6{ Xl -x[) +9(-X n c s )9(\q 1 \ - \X n c i \)9(x 1 - x[) 
-0(-A„cj)0(|A„ a | - \q 1 \)6(x' 1 -x 1 )\$(x,\ n )y{x',\ n )}+h.c., 

where hermitian conjugation is understood in the sense of ©, so that 

P(x,x';q) = iS(x 2 - x' 2 )sgnq 1 6(q 1 (x 1 - x[)) 

N 

^2{tJt>(x,\ n )V(x',\ n ) + t n ^(x,X n )^(x',X n )\ 



N 

X 

, : - J. 



- iS(x 2 - x' 2 )sgnq! 6»(|A nQ | - 

n=l 

x{i„%iA na )<J>(x,A„)$(a;', A,,) + i»0(-giA„cj)$(a;,A»)*(a/,A„)}. 
Again, thanks to ((129(1 the first term cancels out and we get 

N 

P(x,x';q) = ~i6(x 2 - x' 2 ) sgngi ^ 0(|A nQ | - \q x \) 

n=l 

x{t n 0(giA„ a )$(a;,A n )*(a;',A n ) +t„0(-giA„a)$(a;,A n )*(a;', A„)}. (130) 

This suggests the introduction of qi-depcndent solutions of the Nonstationary 
Schrodinger equation 



n (x, qi ) = 0(|A n cj| - |g 1 |){%iA n3 )$(a: ! A„)+f(-9iA„ Q )$(x,A n )}, (131) 
„(*, 9i ) = 0(|A„ a | - |gi|){e(giA„ a )$(a;, A„) + 0(-<ZiA„ 9 )$ (as, A„)}, (132) 



which we call auxiliary Jost solutions. 
Then we rewrite (|13UH as 



jV 



P(x,x';q) = -isgnq 1 6(x 2 - a^) #n(gi)$n(a;, qi)^ n (x', qi) (133) 

n=l 

where 

tfn(ffi) - %iA„ 9 )t» + f(-giA„ 3 )t„. (134) 

Thus the operator P results to be a sum of operators that are different from zero in 
the intervals [— |A„Cs|, |A n cj|] and continuous at x\ — x[. In particular, 

P(q) = for \q x \ > max |A„q|. 

n 

Finally, inserting into 1)68(1 equation (ffjfjfl and equations 182(1 , ((83(1 relating the dressing 
operators to the Jost solutions and the expression obtained for P in 1(133(1 we get the 
completeness relation for the Jost solutions 

— f da <fr(x; a + iq\)t(pt + iqi)^(x'; a + iq\) 
2n J 

N 

- isgngi 2J ^„(<?i)$„(x, gi)*„(a;',gi) = S(xi - x[), for x 2 = x' 2 . (135) 
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4-2. Properties of auxiliary Jost solutions 

Functions & n (x,qi) and ^ n (x,qi) by definitions l|131|) . (|132[) obey the nonstationary 
Schrodinger equation and, correspondingly its dual. Their conjugation property 

$„(x, qi) = Sb n (x,—qi), 

follows from (|86[) . These functions are different from zero on the interval \qi \ < A„q| 
only, that guaranties boundedness of e~ qiXl qi) and e qiXll $> n (x, qi) when x\ goes 

to infinity. As the result the kernel P(x, x '; q) defined by i|31|) is tempered distribution 
with respect to all its variables that proves that operator P(q) belongs to the space of 
operators under consideration. The naturally appeared above piecewise dependence 
of functions $ rl (x, q\) and '^ n {x,qi) on q\ remembers the same, but nontrivial, 
dependence discovered in studying the perturbation of the one- line potential in |18II19| . 
It is necessary to mention that functions <fr n (x,qi) and ^ n (x,qi) are discontinuous at 
qi = 0, while thanks to Ijl29(l P(q) is continuous at this point. 

To study the discontinuity of these functions at qi = we use the standard 
notations (cf. 12411 ) for the right and left limit at q\ = and we write 

$±(x) = 6(±\ n Q)$(x,\ n ) + 9(TKz)$(x,K), 

= 0(±A„ 9 )§O, A„) + 0(TA„3)*(ar, A„), . 

To find relations between these limiting values we have to partially invert relations 
in 11271) and l|128|) . Taking into account that the whole our construction here is 
explicitly symmetric with respect to the A„'s, we can renumber them in a way that 
the first ones for n = 1, . . . , s — 1 are in the upper half plane and the last ones for 
n = s, . . . , N are in the lower half plane. Correspondingly the matrix C decomposes 
as follows 

C=( C Z + Z c _), (136) 

where C± are the matrices, positive and negative definite, introduced in l|117|) . and 
where Z, Z\ which are mutually adjoint, are rectangular matrices just filling the C 
matrix. By introducing the vectors 

$+(A) - . • • , *+ x ) MA) = (M ■ •• , ®- N ) 

d>+(A) fl', ti) *-(A) = (*+, • •• 

and the diagonal matrices 

r = diag{r i }, t = diag{t,}, 

decomposed in the diagonal matrices t± and t± in analogy to ^ 136(1 . equation l|127|l 
can be rewritten as 

*+(*)*+ = i<S> + (X)t + C+tI + ^.(Ajr-ZVj, 

®-{\)t- = i$ + (A)r + ^rl +i$_(A)r_C_rI. 

Solving the last equality for $-(A) and inserting the result into the first one, we get 

N 
n=l 



C" 
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where the matrix C" = Ht^wJ]^ n=1 equals 

(t- 1 y T+ \(C+-ZCZ 1 Z^ 1ZCZ 1 \ ( tUz 1 

(rr^t ) { -iCZ X Z^ —G'Z 1 ){() rZ 

and is one of the limiting values of Ijl34jl : 

4± = 6{±\ n Q)t n + 6(T\n*)i n , 

that in analogy with (|1(J1[I , (|1(J2|) can also be defined as 
d± = res t(k). 

k=A„ !K =Fi| A„ 3 

Thanks to conditions on matrix C given in section l3.4l the matrix C" is hermitian and 
positive. The first property is obvious. In order to check the second one we write for 
an arbitrary vector (w + , i>_ ) 

{v \ v u( c + -zcz 1 z^ lZ CZ l \ ( v + 

= v\C+v+ - [v- -iZ i v + ] f CZ 1 [v-+iZ^v + ], 

that proves that conditions 1)11711 are equivalent to condition C > 0. 

Now relation for the functions ^^(x) follows by complex conjugation, taking into 
account that 



Let us mention that the only dependence on the signs of A„q. enters in the functions 
r(k) in l|98(l . and then in the spectral data (|76JI . All other objects of our constructions 
depend on |A„cj| only. This means that in the case of the zero background potential 
u(x), when both f^ 1 and / ± are identically equal to zero, without loss of generality 
we can chose in addition to IpFfjl A„q > (n = 1, 2, . . . , N). 

On the other side in the case of the nonzero background potential solutions 
corresponding to A„ with opposite sign are different, as was shown in |2()j . 

Finally, let us mention that relation P = P gives the scalar product 

dxi W m {x, qi)® n (x, qi) = — -— - — r — , (137) 

that takes a more complicated form if we turn back to the original solutions. 



4-3. Operators La and M A 

In order to find La we can use the first equality in (|9"U|) applying L to P. Then in 
terms of the hat-kernels we use that &m( x ) 1S annulated by £, so 

N 

L A (x,x';q) ^sgnq 1 S'(x 2 - x' 2 ) ^ tf„(<?i)$ n (x, gi)*„(x', gi). (138) 

n=l 

Now it is easy to check directly that 

PL A = L A P = La j 

as it also follows from from (|9U[) and l|55|) . 

In Sec. 13.31 we proved that in order to obtain Ma it is enough to find, first, a 
general (belonging to our space of operators) solution X of the system l|95fl and then 
to find a special X satisfying 
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The expression for in (|138[) suggests that the general solution X of is 
given by 

N 

)g m (x 2 , x' 2 ; <j)$„(x, qi)^ n (x 1 , q\ ) , (139) 

n=l 

with g m (x 2 ,x' 2 ;q) such that X(x,x';q) 6 S' , but otherwise arbitrary function of the 
written arguments. This can be verified directly by using the representation of P and 
the orthogonality (|137[) . 

Thus we take Ma in the form (|139[) and, according to the discussion at the end 
of Sec. 13.11 we fix the unknown function by requiring that equations l|96|) for X = Ma 
are satisfied. We get 

N 

M A (x,x';q) = - sgn<? x ^ & n (qi){0(x 2 - x' 2 ) + g m (q)}®n(x, gi)$„(x', q x ), (140) 

with g m {q) to be chosen in such a way that Ma(i, x'; q) is bounded at space infinity. 
It results that 

M A (x, x'\ q) = - sgnqi sgn(a; 2 - x' 2 )e~ q ( x ~ x '^ 

N 

x &n(qi)0((q 2 - 2X nmi )(x 2 ~ x' 2 ))$ n (x, qi)^„(x', qi ) (141) 

n=l 

obtained for 

9m{q) = ~0(-q 2 + 2X mSt qi) (142) 

is a well defined bounded operator. 

Since the <£>'s and "f's are bounded we need only to consider the region of variables 
where the exponent e~ q ( x ~ x ) is growing. Since, thanks to the choice l|142|) in each 
term in the r.h.s. of (|14(JH 

-q2(x2 - x' 2 ) < -2\ n ®qi(x 2 - x' 2 ) 
the behaviour of e~ q ( x ~ x ' for qi > cannot be worse than e~ 9l (( ;ri ~ ;l: i) +2A ' l!R ( ;l:2 ~ ;r 2)) 
and, then, since \q\\ < |A„q|, no worse than 9(±\ n c ! )e Tes< ' :Xn ' >< ' x ~ x '. But in the interval 
< qi < |A„q| where <£>„ and ^> n are different from zero 

9(±X nQ )e^ x ^ x - x '^ n (x, qi)* n (x', Ql ) 

= 6(±\ n , 3 )e~ u ^ x ^ x - x '^x(xAn ;il +i\Xn*\)Z{x', A„ K +i\\ n z\) 

is bounded. Analogously for qi < 0. This proves that Ma is bounded. 

We are left with the explicit construction of (MC • Form (|70|l and the conjugation 
properties of the Jost solutions we have 

CMC 1 = vtujMvt ] u) 

and, then, from lfT7jl. lf2U|) and [Rjfij) 

CMC} = VTMqu. 
In the p-space this reads 

(CMC%;q)= f dp'u(p-p';q 1 +p' 1 ) J^ 1+ ( P ' l \ ^ (p'iqi) 

J q2 +p 2 - (qi +Pi) 
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and shifting p', that is naming p' + qsR = a, 

((M^)(p:q) = [ dav{p- a + c^;ai+iqi) ^(gl+jgij _ _ ^ 

Recalling (|37|l and inverting ijBJ we derive for C-^C m the x-space 

(C^C t ){x 1 x';q) = ——- / da — . x(x; ai + %)£(x'; ai + %), 

(2ny J a 2 + tq 2 - {ai + tqi) 2 

and, finally, integrating over ct<2, and summing up Ma as indicated in (|141|1 we get for 
the hat version of the resolvent M 



x — IdaiO {{q 2 - 2a 1 q 1 ){x 2 - x' 2 )) t{a\ + %)$(x; ai + ai + %) 



AT 



M(x, x';q) = sgn(x2 - x' 2 ) 
1 

2nij 

- sgngi sgn(x 2 - x' 2 ) ^ i?n(gi)^((£?2 - 2A„sRgi)(:E2 - x 2 ))$ n (x, gi)$„(x', 91). (143) 
n=l 

It is worthwhile to note that in the r.h.s. of this formula for M the poles of £(qi) at 
qi = A„ and at qi = A„ in the first term generate discontinuities at q\ = ±A n a which 
cancel exactly the discontinuities of the second term, due to the fact that the auxiliary 
Jost solutions $„(q , i), * n (gi) are identically zero outside the strip \q\\ < |A n a|. 

By using the reduction procedures indicated in 1-iL'l) and (|33|l from this explicit 
expression for the resolvent one can derive the generalization of the standard 
Green's function for Jost and advanced/rearted solutions to the case of N solitons 
superimposed to a generic background. Thanks to the fact that they are derived from 
the resolvent they result both to be bilinear in terms of the Jost and auxiliary Jost 
solutions. The use of these Green's functions and other ones suggested by studying 
the properties of the resolvent is crucial in extending the 1ST to potentials obtained by 
perturbing the potential considered in this paper by adding to it an arbitrary smooth 
rapidly decaying function of both spatial variables. This is performed in a forthcoming 
paper following the method developed in |18j . 
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